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ABSTRACT 

IWalter et al.l ()201 2^ have recently identified the J=6-5, 5-4, and 2-1 CO rotational 
emission lines, and [Cii] fine-structure emission line from the star-forming interstellar 
medium in the high-redshift submillimeter source HDF 850.1, at z = 5.183. We employ 
large velocity gradient (LVG) modeling to analyze the spectra of this source assuming 
the [Cii] and CO emissions originate from (i) separate unvirialized regions, (ii) sep- 
arate virialized regions, (in) uniformly mixed unvirialized region, and (iv) uniformly 
mixed virialized regions. We present the best fit set of parameters, including for each 
case the ratio a between the total hydrogen/helium gas mass and the CO(l-O) line lu- 
minosity. We also present computations of the ratio of H2 mass to [Cii] line-luminosity 
for optically thin conditions, for a range of gas temperatures and densities, for direct 
conversion of [Cii] line-luminosities to "dark-H2" masses. For HDF 850.1 we find that 
a model in which the CO and C"*" are uniformly mixed in gas that is shielded from UV 
radiation, requires a cosmic-ray or X-ray ionization rate of ^ ~ 10"^"^ s~^, plausibly 
consistent with the large star- formation rate (~ 10'^ yr~^) observed in this source. 
Enforcing the cosmological constraint posed by the abundance of dark matter halos in 
the standard ACDM cosmology and taking into account other possible contributions 
to the total gas mass, we find that three of these four models are less likely at the 
2(T level. We conclude that modeling HDF 850. 1's ISM as a collection of unvirialized 
molecular clouds with distinct CO and C+ layers, for which a = 0.6 M0 (K km s^^ 
pc^)~^ for the CO to H2 mass-to-luminosity ratio, (similar to the standard ULIRG 
value), is most consistent with the ACDM cosmology. 

Key words: cosmology: theory - galaxies: high-redshift - galaxies: ISM - galaxies: 
mass function. 



1 INTRODUCTION 

Observations of high-redshift CO spectral line emissions 
have greatly increased our knowledge of galaxy assembly 
in the early Universe. At redshifts z~2, this includes the 
discovery of turbulent star-forming disks with cold-gas mass 
fractions and star-for matio rates significantly larger than in 
present day galaxies ( Daddi et al. 2010l: Genzel et al]|2012l : 
ikagnelli etal] I2OI2I : iTacconi et al.ll20ld . |2012| ). The high 
star-formation rates are correlat ed with large gas masses 
and luminous CO emission lines (|Kennicutt fc Evansll20l3 ) 
observable to very high redshifts. A prominent example is 
the luminous submillimeter and Hubble-Deep-Field source 
HDF 850.1, which is at a redshift of z = 5.183 as determined 
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by recent detections of CO(6-5), CO(5-4), and CO(2-l) ro- 
tational line emissi ons, and also [Cii] fine-structure emis- 
sion in this source ([Walter et al.ll2012l '). The high redshift 
of HDF 850.1 offers the opportunity of setting cosmological 
constraints on the conversion factor from CO line luminosi- 
ties to gas masses, via the implied dark matter masses and 
the expected cosmic volume density of halos of a given mass. 
Such analysis is the subject of our paper. 

CO emitting molecular clouds, which provide the raw 
material for star formation, are usually assumed to have 
undergone complete conversion from atomic to molecular 
hydrogen. However, since H2 has strongly forbidden rota- 
tional transitions and requires high temperatures (~500 K) 
to excite its rotational lines, it is a poor tracer of cold (< 
100 K) molecular gas. Determining H2 gas masses in the 
interstellar medium (ISM) of galaxies has therefore relied 
on tracer molecules. In particular, ^^CO is the most com- 
monly employed tracer of ISM clouds; aside from being the 
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most abundant molecule after H2, CO has a weak dipole 
moment (^e = 0.11 Debye) and its rotational levels are thus 
excited and thermalized by collisi ons with H2 at relatively 
low m olecular hydrogen densities (jSolomon fc Vanden Bou3 
I2OO5I ). 

The molecular hydrogen gas mass is often obtained from 
the CO luminosity by adopting a mass-to-luminosity conver- 
sion factor a = MhiI1^'(_ 



CO(l-O) 



between the H2 mass and 



the J =1-0 115 GHz CO rotational transition (jSolatto et alJ 
I2OI3I ). The value for a has been empirically calibrated for 
the Milky Way Galaxy by three independent techniques: (i) 
correlation of optical ext inction with CO column densities in 
interstellar dark clouds (|DickmarJll978l ): (ii) correlation of 
gamma-ray flux with the CO line flux in the Galac tic molec- 
ular ring (Bloem en et al. l ll986l : IStrong et al.|[l98^ ): and (m) 
observed relations between th e virial mass and CO line lu- 
minosity for Galactic GMCs (|Solomon et al.lll987h . These 
methods have all arrived at the conclusion that the con- 
version factor in our Galaxy is fairly constant. The stan- 
dard Galactic value is a = 4.6 M0/(K km"^ pc^). Subse- 
quent studies of CO emission from unvirialized regions in 
Ultra-Luminous Infrared Galaxies (ULIRGs) found a sig- 
nificantly s maller ratio. For such syste ms, a = 0.8 M0/(K 
km~^ pc^) (|Downes fc Solo mon 1998). These val ues have 
been adopted by many |Sandcrs et al. 1988 : Tinnev et al.l 
I1990I : IWang et al.lll99]J : IWaher et alTl2012l ) to convert CO 
J=l-0 line observations to total molecular gas masses, but 
without consideration of the dependence of a on the av- 
erage molecular gas conditions found in the sources being 
considered. Since all current observational studies of a leave 
its range and dependence on the average density, temper- 
ature, and kinetic state of the molecular gas still largely 
unexplored, its applicability t o other systems in the loca l or 
distant Universe is uncertain (jPapadopoulos et al.ll20l3 ). 

In this paper, we estimate the CO emitting gas masses 
in HDF 850.1 using the large- velocity-gradient (LVG) for- 
malism to fit the observed emission line spectral energy dis- 
tribution (SED) for a variety of model configurations, and 
for a comparison to the Galactic and ULIRG conversions. 
In §2 we outline the details of the LVG approach, including 
an overview of the escape probability method and a deriva- 
tion of the gas mass from the line intensity of the modeled 
source. In §3, we present the best fit set of parameters that 
reproduce HDF 850. 1's detected lines and calculate the cor- 
responding molecular gas mass, assuming the CO and [Cii] 
emission lines originate from (i) separate unvirialized re- 
gions, (ii) separate virialized regions, ('m^uniformly mixed 
unvirialized regions, and (iv) uniformly mixed virialized re- 
gions. The inferred gas masses enable us to set lower-limits 
on the dark-matter halo mass for HDF 850.1. In §4 we com- 
pare the estimated halo-masses to the number of such ob- 
jects expected at high redshift in ACDM cosmology, and 
show that models with lower values of a are favored. We 
conclude with a discussion of our findings and their impli- 
cations in 55. 



2 LARGE VELOCITY GRADIENT MODEL 

We start by describing our procedure for quantitatively an- 
alyzing the [Cii] and CO emission lines detected at the po- 
sition of HDF 850.1 using the large velocity gradient (LVG) 



approximation. We consider a multi-level system with pop- 
ulation densities of the ith level given by rii . The equations 
of statistical equilibrium can then be written as : 



(1) 



where I is the total number of levels included; since the set 
of I statistical equations is not independent, one equation 
may be replaced by the conservation equation 



ntot 



(2) 



where ntot is the number density of the given species in all 
levels. In our application, ntot = n co- Following the no- 
tation of IPoelman fc Spaanj l|2005l ). Rtj is given in terms 
of the Einstein coefficients, Aij and Bij, and the coUisional 
excitation (i < j)and de-excitation rates (i > j) dj: 



Rii 



A, j+B,j{J^,)+C,j, {i>j) 

B, ,{J^j)+a,, ii<j) 



(3) 



where (Jij) is the mean intensity corresponding to the tran- 
sition from level i to j averaged over the local line profile 
function (f){h'ij). The total coUisional rates dj depend on 
the individual temperature-dependent rate coefficients and 
collision partners, usually H2 and H for CO 

The difficulty in solving this problem is that the mean 
intensity at any location in the source is a function of the 
emission and varying excitation state of the gas all over the 
rest of the source, and is thus a nonlocal quantity. To obtain 
a general solution of the coupled sets of equations describ- 
ing radiative transfer and statistica l equi librium, we adopt 
the approach devel oped bylSobplevI (|l96d ) and extended by 
ICastoil l|l970t ) and ILucvI l|l97ll ) and assume the existence 



of a large velocity gradient in dense clouds. This assump- 
tion is justified given the interstellar molecular line widths 
which range from a few up to a few tens of kilometers per 
second, far in excess of plausibl e thermal velocities in the 
clouds (|Goldreich fc Kwanlll973 ). They suggest that these 
observed velocity differences arise from large-scale, system- 
atic, velocity gradients across the cloud, a hypothesis that 
lies in accord with the constraints provided by observation 
and theory. 

In the limit that the thermal velocity in the cloud is 
much smaller than the velocity gradient, the value of {Jij) 
at any point in the cloud, when integrated over the line pro- 
file, depends only upon the local value of the source function 
and upon the probability that a photon emitted at that point 
will escape from the cloud without further interaction. Thus 
(Jij) becomes a purely local quantity, given by: 



{,hj) = {l-P^j)S.j+P.,B{u,j,TB) 



where Sij is the line source function, 
_ 2hu% ^ g,nj 



1)- 



(4) 



(5) 



assumed constant through the medium. In this expres- 
sion Qi and Qj are the statistical weights of levels i and 
j respectively, Pij is the "photon escape probability", and 
Bij (vij , Tb ) is the background radiation with temperature 
Tb- In our models we set Tb to the CMB temperature of 
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Table 1. CO and [Cll] Line Observations of HDF 850.1 



Line 




Integrated Flux Density 


Luminosity 


Intensity 




[GHz] 


[Jy km s^-"-] 


[10^0 Kkms-i pc2] 


[erg(s cm^ str)"-"-] 


CO{2-l) 


37.286 


0.17±0.04 


4.1±0.9 


(2.16±0.51)xl0-^ 


CO(5-4) 


93.202 


0.5±0.1 


1.9±0.4 


(1.59±0.32)xl0-« 


CO(6-5) 


111.835 


0.39±0.1 


1.0±0.3 


(1.49±0.82)xl0-* 


[Cii] 


307.383 


14.6±0.3 


5.0±0.1 


(1.5±0.03)xl0-s 



Detection of four lines tracing the star-forming interstellar medium in HDF 850.1 l lWalter et al.ll2012l'l 



16.9 K at z = 5.183, an d we ignore contributions from warm 
dust (|da Cunha et all 12013). 

For a spherical homogenous collapsing cloud, the prob- 
ability that a photon emitted in the transition from level i 
to level j escapes the cloud is given by 



(6) 



A. = ^ 

where Tij is the optical depth in the line. 



^ _ Ajj c ntot njQi QiTii 
'■^ Stt vf- dv/dr ntotOj Qirij 



(7) 



The equations of statistical equilibrium are therefore re- 
duced to the simplified form 

I I I 

^(jiiCij — UjCji) + Tii ^ AijPij — ^ UjAjiPji = (8) 



and can be solved through an iterative process to give the 
fractional level populations ru/ntot (for a given choice of den- 
sities for the collision partners, usually H2 but also H, and ki- 
netic temperature Ttin). Assuming the telescope beam con- 
tains a large nu mber of these identical ho mogeneous collaps- 
ing clouds fe.g.. lHailev-Dunsheathir2009l ). the corresponding 
emergent intensity of an emission line integrated along a line 
of sight is then simply 



j-ii — — : ^ijPij^ytot 



47r ntot' 



^^A,,A:.X^ (9) 
47r n 



where x is the abundance ratio, x = ntot I n, where n is the 
total hydrogen gas volume density (cm~^) and N is the hy- 
drogen column density (cm~^. 

Given a series of observed lines with frequency Vij , one 
can identify a set of characterizing parameters that best re- 
produces the observed line ratios and intensity magnitudes; 
among these parameters are the cloud's kinetic temperature 
(Tkin), velocity gradient (dv/dr), gas density n and colli- 
sion partner (H and H2) gas fractions, abundance ratio (x), 
and column density (A'^). For a spherical geometry, this col- 
umn density can then be further related to the molecular 
gas mass of the cloud in the following way 



M^-n-R^fimN' 



(10) 



where the factor p = 1.36 takes into account the helium 
contribution to the molecular weight and 7? = Da9/2 is the 
effective radius of the cloud, with Da being the angular di- 
ameter distance to the source and 6 the beam size of the 
line observations. is defined in terms of the column den- 
sity obtained from the LVG calculation, A''' = A'^(l + z)* 
where the (1 -I- z)* multiplicative factor reflects the decrease 



in surface brightness of a source at redshift z in an expand- 
ing universe. 

In the case where the emitting molecular clouds are 
gravitationally bound, applying the virial theorem to a ho- 
mogenous spheri cal body yields the following constraint 
(|Goldsmithll200ll ). 



dv / G-K^mn 
dr ~ V 15 



(cm~3) 



km s ^ pc ^ 



(11) 



where a = 7.77x10"^ if n = and a = 5.50x10^^ if n 
— riH ■ We note that the velocity gradient is inversely pro- 
portional to the dynamical time scale. Therefore, in models 
where the clouds are assumed to be virialized, dv/dr and n 
are no longer independent input parameters of the model, 
but rather vary according to equation (11). 

For the optically thick ^^CO J=l-0 transition line (/3 ~ 
1/r), carrying out the LVG calculations with this additional 
virialization condition leads to a simple relation between the 
gas mass and the CO (1-0) line luminosity, 



M. 



OLCO = 



H2 



= 8.6 



1 (c 



-^'cO(l-O) T(!a:c/(K) 



-M0(K km s ^ pc 



(12) 

where the excitation temperature T^^c ~ T^in when the 
emission line is thermalized. (In this expression we assume 
complete conversion to H2 so that the gas mass is the 
H2 mass.) Empirically, the Galactic value of this mass-to- 
luminosity ratio for virialized objects bound by gravitational 
forces is g = Mg,/L'(7oa -oi = 4.6 Mq (K km s"^ pc^)"^ 
(|Solomon fc Barrett|[l99ll ). 



3 ANALYSIS OF HDF 850.1 
3.1 Properties of the Galaxy 

The detection of HDF 850.1, the brightest submillimetre 
source found in the Hubble Deep Field at a wavelength of 
850 micrometers, was confirmed bv lHughes et al. 1998. Since 
then, a full-frequency scan towards the Hubble Deep Field 
(beam size of 2.3", n « 10"^° radians) using the IRAM (In- 
stitut de Radioastronomie Millimetrique) Pleateau de Bure 
Interferometer and the National Radio Astronomy Observa- 
tory (NRAO) Jansky Very Large Array has detected three 
CO lines, identified as the CO(2-l) 230.5 GHz, CO(5-4) 
576.4 GHz, and CO(6-5) 691.5 GHz rotational transitions. 
[Cii] 1900.1 GHz (158 /im), one of the main cooling lines of 
the star-forming in stellar medium, has also been detected 
(Table 1). These lines have placed HDF 850.1 in a galaxy 
over density at z = 5.183, a redshift higher than those of 
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Tkin [K] 

Figure 1. The gas mass to [Cll] luminosity ratio for a C"*" to H2 abundance ratio of ^^+=10"'* as a function of the kinetic temperature 
Tfei„ at a molecular hydrogen number density of rijj^ = 10, 10^, 10'', and 10* cm~^. Calculations were made in the optically thin 
regime, assuming the fine-structure transition J=3/2— >l/2 is due solely to spontaneous emission processes and collisions with ortho- and 
para-H2 . 



most of the hu ndreds of submillim etre-briRht galaxies iden- 
tifie d thus far (IWalter et al.l I2OI2I ) . 

IWalter et alj (|2012l )" used an LVG model to character- 
ize the CO spectral energy distribution of HDF 850.1 and 
found that the observed CO line intensities could be fit with 
a molecular hydrogen density of 10^'^ cm~^, a velocity gra- 
dient of 1.2 km s^^pc"^, and a kinetic temperature of 45 
K. Then, assuming a = 0.8 M0(K km s~^pc'^)~^ as for 
ULIRGs, they used the 1-0 line luminosity inferred from 
their LVG c omputation to infer that M g, = 3.5x10^° Mq. 
However, as IPapadopoulos et"aLl (|2012l 'l argue, adopting a 
uniform value of a for ULIRGs neglect its dependence on 
the density, temperature, and kinematic state of the gas; 
this may Umit the applicability of computed conversion fac- 
tors to other systems in the local or distant Universe. 

Here, we bro aden the LVG analysis carried out in 
IWalter et aLll2012l and use the LVG-modeled column den- 
sity to estimate the total gas mass of the source. We present 
several alternative models, each subject to a slightly differ- 
ent constraint. In particular, we first consider the case where 
the CO and [Cii] lines originate from different regions of the 
molecular cloud and then consider models in which the CO 
molecules and C"*" ions are uniformly mixed, such that the 
line emissions originate in gas at the same temperature and 
density. In both instances, we perform our LVG computa- 
tions assuming (a) gravitationally-unbound clouds and (b) 
virialized clouds for which the virialization condition (equa- 
tion (11)) has been imposed. 

To carry out these computatio ns, we use the LVG ra- 
diative transfer code described in iDavies et af] I2OI2I . En- 
ergy levels, line frequencies and Einstein A coefficients 



are taken from the Cologne Database for Molecular Spec- 
troscopy (CDMS).The excitation and deexcitation rates of 
the CO rotational l evels that ar e indu ced by collisions with 
H2 are taken from lYang et al.l (|2010|) while the C"*" coUi- 
sion al rate coefficients come from Flower fc Launavl (|l977l ) 
and iLaunav fc Rouefj (|l977l ). 



3.2 Separate CO, C+ Unvirialized Regions 

We first consider a model in which the CO and [Cii] emission 
lines detected at the position of HDF 850.1 originate in sep- 
arate regions of the molecular gas cloud, regions which are 
not necessarily at the same temperature and number den- 
sity. This picture is consistent with the standard structure of 
PDRs in which there is a layer of almost totally ionized car- 
bon at the outer edge, intermediate regions where the carbon 
is atomi c, and internal regions where the carbon is locked 
into CO (ISternberg fc Dalgarnd[l995l : iTielens fc Hollenbaciil 
ll985l : IWolfire et allboiol v" In this picture then, the hydrogen 
is fully molecular in the CO emitting regions. 

We apply our LVG model to the observed CO line inten- 
sities, without the additional virialization constraint given 
by equation (11). We find that the CO SED is best fit with a 
molecular hydrogen density of lO^ '^ cm~^, a kinetic temper- 
ature of 85 K, a velocity gradient dv/dr = 5 km s~^pc~^. 



and a CO to H2 abundance ratio of X.qq = 10~ . The col- 
umn density that yields the correct line intensity magnitudes 
is 5.1x10^* cm~^, corresponding to a molecular hydrogen 
gas mass of ~ 2.54x10^° Mq. This est imate of the 

gas m ass is nearly 30% less than obtained b\ i Walter et al.l 
(I2OI2I) by applying the H2 mass-to-CO luminosity relation 
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with the typically adopted conversion factor for ULIRGs, a 
= 0.8Mq(K km s-^ pc2)-\ Given our inferred H2 eas mass 
and predicted CO(l-O) line luminosity from our LVG fit, we 
find that a = 0.6 M0(K km s~^ pc^)~^, in this model. 



3.3 Separate CO, C+ Virialized Regions 

For self-gravitating clouds in virial equilibrium, the veloc- 
ity gradient is no longer an independent input parameter of 
the LVG model, but varies with according to equation 
(11). We calculate a three-dimensional grid of model CO line 
SEDs, varying Tkin, , and xco over a large volume of the 
parameter space. We find, under this virialization constraint, 
that the observed CO lines are fit best with a temperature of 
100 K, an abundance ratio of xco ~ 10~*'^, and a molecular 
hydrogen number density of 10^'^ cm~'^ (with a correspond- 
ing velocity gradient of 0.20 km s~^pc~^). For this set of 
parameters, the beam-averaged H2 column density is well 
constrained to be N_r/2 « 5.2x10^® cm~'^, giving an associ- 
ated molecular gas mass of Mh2 ~ 2.59x10" M0. The H2 
mass to CO luminosity conversion factor obtained in this 
model is a = 5.9 M0(K km s"'^ pc'^)"^, a value similar to 
the Galactic conversion factor observed for virialized molec- 
ular clouds in the Milky Way, suggesting that HDF 850.1 
may have some properties in common with our Galaxy. The 
factor-of-ten increase in the inferred mass when imposing 
the virialization condition is a measure of the uncertainty in 
the CO to H2 mass conversion even for a standard "FDR 
picture" . 



L[C7i] = 1.1x10^" 



Lq, is found to be « 5.2x10^° 



3.4 Optically thin [Cii] 

In the two models above, where the CO and [Cii] lines are as- 
sumed to be emitted from separate regions of the molecular 
clouds, the single detected ionized carbon line is insufficient 
in constraining the parameters of the LVG modeled [Cii] re- 
gion. We thus consider the optically thin regime of the [Cii 
line (/3 — 1), such that 



Mh2 _ fJ-mH2{^ + zY 



L 



(13) 

where Xc+ ^'^ ^2 abundance ratio (fixed at a value 

of 10~* in these calculations) and Xi represents the frac- 
tion of ionized carbon molecules in the ith energy level. As- 
suming that the fine-structure transition J=3/2— >-l/2 is due 
solely to spontaneous emission processes and collisions with 
ortho- and para-H2 (radiative trapping can be neglected in 
the optically thin regime), the equations of statistical equi- 
librium reduce to a simplified form and can be solved to 
obtain 2:(j^3/2) as a function of temperature and H2 num- 
ber density. In Fig. 1, we have plotted the resulting mass- 
to-luminosity ratio, MH2/ L[cii], as a function of the kinetic 
temperature for several different values of nH2- Given the 
high [Cii] /far- infrared lumino sity ratio of Lfcin /LpiR = 
(1.7±0.5)xl0~^ in HDF 850.1 (|Waher et al.ll2012l i. it is rea- 
sonable to assume the ionized carbon is emitting efficiently 
and to thus consider the high temperature, large number 
density limit {Tkin — 500 K, n^r^ ~ 10* cm~^). In this limit, 
the mass-to- luminosity ratio is ~ 1.035 M0(K km~^ pc^)~^ 
and the corresponding molecular gas mass of the C^ region, 
using the detected line intensity of the ionized carbon line 



3.5 Uniformly Mixed CO, C"*" Unvirialized Region 

We next consider models in which the CO molecules and C^ 
ions are mixed uniformly, such that the corresponding line 
emissions originate in gas at the same temperature, den- 
sity, and velocity gradient. These models resemble condi- 
tions found in UV-opaque cosmic-ray dominated dark cores 
of interstellar molecular clouds where the chemistry is driven 
entirely by cosmic-ray ionization. For these conditions, the 
density fractions rii/n for CO and C"*" depend on a single pa- 
rameter, the ra tio of the c loud density to the cosmic-ray ion- 
ization rate (|Boger fc Sternbera.2005. ). For Galactic con- 
ditions with ~10~^^ s~^, the C^ to CO ratio is generally 
very small. However, in objects such as HDF 850.1 where 
the ionization rate may be much larger due to high star for- 
mation rates, this ratio may be enhanced significantly. 

The LVG fit for the observed set of line intensity ra- 
tios, {I[CII]/]^CO{2-l), l\CII]/IcO(5-'l), IiCII]/IcO(6-5)} = 

{7.08±1.67, 0.96±0.19, 1.03±0.2}x 10^ (|Walter et al.ll2012l ). 
yields Tki„ = 240 K, n = 10^'^ cm"^ and a C+ to CO abun- 
dance ratio of 12. To estimate the ionization rate required 
to produce this C^ /CO a bundance ratio, we employed the 
iBoger fc Sternberg! (|2005l ) chemical code that captures the 
C^ - C - CO interconversion in a purely ionization-driven 
chemical medium. We solve for the corresponding gas phase 
metallicity and cosmic-ray ionization rate pairs {Z, C^) that 
yield this relatively high C""" to CO ratio for a cloud with 
temperature T=240 K and density nH2=10^ '' cm~^. 

Assuming a Galactic metallicity of nearly one, (a fairly 
reasonable assumption given the high star formation rate 
observed in HDF 850.1), the cosmic-ray ionization rate re- 
quired to achieve this high C^ to CO ratio is of order ^ 
10" 



s~^ (Figure 2). Compared to the cosmic-ray ioniza- 
j tion rate (-^lO"'^® s~^) estimated for the Milky Way, the 
rate found in HDF 850.1 is indeed significantly enhanced. 
This finding is consistent with the fact that HDF 850.1 has 
a star-formation rate of 850 solar masses per year, a value 
which i s larger than the m easured Galactic SFR by a factor 
of 10^ l|Walter et al.ll2012l ). 

Furthermore, for such a high cosmic-ray ionization rate 
of this magnitude, the hydrogen is primarily atomic for the 
implied LVG gas density. We find that the ratio of atomic 
hydrogen to molecular hydrogen in the interstellar clouds 
is n_ff/nH2 ~ 3 (Figure 2). We thus replace H2 with H as 
the dominant collision partner. Given the uncerta in CO-H 
rotational excitation rates fsee IShepler et al]|2007l ). we as- 
sume that the rate coe fficients are equal to the rates fo r col- 
lisions with ortho-H2 (jFlower fc Pineau Des Foij|2010l ). We 
find that the observed CO and [Cii] lines together are best 
fit with a temperature of 240 K, a velocity gradient of 1.5 



-1 
n3.3 



= 10 cm , and an abundance ratio (relative to H) of 
7.4x10"^ and 6x10"^ for C+ and CO respectively. The col- 
umn density that yields the correct line intensity magnitudes 
is 3.3x10^® cm~'^, corresponding to an atomic hydrogen gas 
mass of ^ 8.15x10^° Mq. The molecular hydrogen mass 
is then Mh^ = 2{nH2/nH)MH ~ 4.99x10^° M0 and the to- 
tal gas mass estimate is Mgas ~ 1.31x10^^ M0. The corre- 
sponding conversion factor in this model is a = 5.3 M0(K 
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Figure 2. Dependence of the density ratios, Uqj^ jnco (solid line) and nH^jriH (dashed line), on the H2 ionization rate ^ at a fixed 
solar metallicity ^=1, for our best-fit LVG parameters Tfci„ = 240 K and njj^ = lO^ '^ cm~^). As C, increases, the abundance of C+ 
relative to CO in the cosmic-ray dominated dark cores of interstellar clouds grows while that of H2 to atomic hydrogen decreases. The 
desired value nQj^. jnco ~ 12 is obtained for 10~^^ at which point nu Sn^f,, . 



km s pc ) , where a is now defined as the total gas mass 
to CO luminosity ratio. 



3.6 Uniformly Mixed CO, C+ Virialized Region 

In the case where we assume self-gravitating clouds in virial 
equilibrium, we again find a best fit model with a to 
CO ratio of Xc+ /Xco ~ indicating that atomic hydro- 
gen is the dominant collision partner in the LVG calcula- 
tions. We thus calculate a three-dimensional grid of model 
CO and [Cii] fines, varying T^i^, nn, X[ci/]' ^^'^ Xco- 
observed set of CO and [Cii] fines, assumed to have been 
emitted from the same region, are fit best with a tempera- 
ture of 230 K, an atomic hydrogen number density of 10'^ 
cm~'^ (with a corresponding velocity gradient of 0.17 km 
pc~^) and an abundance ratio of 7.3xl0~^ and 7xl0~^ 
for C"*" and CO respectively. For this set of parameters, the 
beam-averaged H column density is well constrained to be 
N_tf ~ 8.5x10^® cm~^. This corresponds to an atomic and 
molecular gas mass of Mh « 2.12x10" Mq and Mh^ ^ 
1.44x 10^^ M0 respectively, yielding a total gas mass esti- 
mate of Mgas ~ 3.56x10" Mq. The ratio of the total gas 
mass to the CO luminosity in this model is a = 8.5 M0(K 
km s"^ pc2)-\ 



4 COSMOLOGICAL CONSTRAINTS 

Our inferred gas masses enable us to set cosmological con- 
straints. For a particular set of cosmological parameters, 
the number density of dark matter halos of a given mass 



can be inferred from the halo mass function. The Sheth- 
Tormen mass function, which is based on an ellipsoidal col- 
lapse model, expresses the comoving number density of halos 
n per logarithm of halo mass M as. 



nsriM) 



dn 



dlogM 



A 



2a pr, 



dv 



■K M dlogM 



(a//2 



(14) 

where a reasonably good fit to simulations can be ob- 
tained by scttmg A = 0.322, a = 0.707, and p = 0.3 
(Shcth & TormeHlHii). Here, pm is the mean mass den- 
sity of the universe and v = Scrit{z) /o-{M) is the number of 
standard deviations away from zero that the critical collapse 
overdensity represents on mass scale M. Integrating this co- 
moving number density over a halo mass range and volume 
element thus yields TV, the total number of halos observed 
within solid angle A with mass greater than some Mh and 
redshift larger than some 2, 



N{z, Mu) = 



Ho 



A / dz- 



' Ml, 



dlogM 



dn 



dlogM 



(15) 



where Ho is the Hubble constant. Da is the angular diam- 
eter distance, and Qm and Qa are the present-day density 
parameters of matter and vacuum, respectively. 

Given the detection of HDF 850.1, we can say that out 
of the hundreds of submillimetre-bright galaxies identified so 
far, at least one has been detected in the Hubble Deep Field 
at a redshift z > 5 with a halo mass greater than or equal 
to the halo mass associated with this source. This obser- 
vation, taken together with the theoretical number density 
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Table 2. LVG Model: Best Fit Parameters and Results 



Model Parameters 


Separate, 


Separate, 


iviixed, 


iVlixea, 




unvirialized 


virialized 


unvirialized 


virialized 




85 


100 


240 


230 


dv/dr [km s"^ pc^^] 


5.0 


0.2 


1.5 


0.17 


[cm-3] 


103.2 


102.8 


102.8 


102.5 


UH [cm"3] 






103.3 


10^ 


Xco/H2 


10-3.8 


10-4.5 


10-4.7 


10-4.7 


Xco/H 






10-5.3 


10-5.2 




0.5 


5.2 


1.0 


2.9 


Nh [IQi^ cm-2] 






3.3 


8.5 


Model Results 


[Mq] 


2.54x10^° 


2.59x10" 


4.99x1010 


1.44x10" 


Mh [Mq] 






8.15x101" 


2.12x10" 


Mgas [Mq] 


2.54x101" 


2.59x10" 


1.31x10" 


3.56x10" 


LcO{i-0) [(K km s-i pc2)] 


3.95x10"' 


4.41x10"' 


2.49xl0i" 


4.18x101" 


[M0(K km s-i pc2)-i] 


0.6 


5.9 


5.3 


8.5 



a here is defined as the ratio of total gas mass to CO(l-O) luminosity, a = Mgas/L'cO(i-0) • the 
models where the CO and [Cll] lines are assumed to be originating from separate regions, Mg^s = M//^ 
since estimates of the atomic hydrogen gas mass could not be obtained via the LVG calculations. In the 
models where the CO molecules and C+ ions are assumed to be uniformly mixed, the total gas mass is the 
sum of the molecular and the atomic gas masses, Mgas = M//2+Mif. 



predicted by the Sheth-Tormen mass function, constrains 
the halo mass associated with HDF 850.1 in the following 
way, 

iV(5, Afh.™„) > 1 (16) 

where the solid angle co vered by the Hub ble Deep Field 
is Ahdf = 4.4 arcmin^ (|Flvnn et al.ll200ll ') and a ACDM 
cosmology is assumed with Hp = 70 km s" ^ Mpc"!, f^A 
= 0.7, and r2„ = 0.3 (|Komatsu et al.ll201ll '). Mh,mm, the 
minimum inferred halo mass for HDF 850.1, is related to 
the halo's minimum baryonic mass component, a quantity 
derived in §3 via the LVG technique, in the following way, 

Mfi,min = -TT—Mb,min (17) 
lit 

where the baryonic and the total matter density parameters 
are Qb = 0.05 and flm ~ 0.27 respectively. Each model's es- 
timate of the minimum baryonic mass associated with HDF 
850.1 corresponds to an estimate of N, the number of ha- 
los expected to be observed at a redshift z > 5 with mass 
Mh.min- Given that one such source has already been ob- 
served, an atomic physics model that yields a number A'' < 
1 can be ruled out at a confidence level P = 1 - N; con- 
versely, a model with a baryonic mass estimate that yields 
N ^ 1 is consistent with the cosmological constraint given 
by equation (15) and remains a viable model candidate. 

The confidence with which models can be ruled out on 
this basis is plotted as a function of the minimum baryonic 
mass estimated by the model (solid curve in Figure 3). To 
check the consistency of the four LVG models considered in 
§3 with these results, dashed lines representing the masses 
derived from each model are included in the plot (upper left 
panel). Assuming the CO and [Cii] molecules are uniformly 
mixed in virialized clouds results in a baryonic mass Mb^min 
= 3.56x10" Mq. Such a source, with a corresponding halo 
mass Mh ^ 1.9x10^^^ Mq, is expected to be observed 4% of 



the time; this model can thus be ruled out at the 2a level 
(dotted) . The models which postulate separate virialized re- 
gions (dash-dot) and mixed unvirialized regions (solid) can 
be ruled out with relatively less certainty, at the 1.6ct and 
0.5a levels, respectively. On the other hand, modeling the 
CO and [Cii] emission lines as originating fr om separate, 
unviri alized regions (dashed), as was done bv lWalter et al.l 
(|2012l ), results in a minimum baryonic mass which is con- 
sistent with the constraint posed by equation (15). In this 
case, Mb,min ~ 2.54x10^" Mq, and we expect to find A'' > 
1 halos with a halo mass ^ 5.4x10^^ M0 at a redshift z 
^ 5. 

These baryonic masses, obtained via the LVG method, 
represent conservative estimates of the total baryonic con- 
tent associated with HDF 850.1. In particular, mass contri- 
butions from any ionized gas or stars residing in the galaxy 
are not taken into account, and in the models where a 
layered cloud structure is assumed, the atomic hydrogen 
mass component is left undetermined. We therefore con- 
sider the effects on the predicted number of observed halos, 
N(Mh,min), if the minimum baryonic mass derived for each 
model is doubled. Using these more accurate estimates of 
HDF 850. 1's baryonic mass, we find that the two models 
in which the virialization condition is enforced can be ruled 
out at the > 2a level (upper right panel). 

For comparison, we also consider the gas masses implied 
by the CO (1-0) line luminosities predicted by the two mod- 
els where distinct CO and C^ layers are assumed (lower left 
panel). Adopting a CO-to-H2 conversion factor of a = 0.8 (K 
km s~i pc^)"^, we find that enforcing the cosmological con- 
straint posed by the abundance of dark matter halos does 
not rule out either model, even if the minimum baryonic 
masses are doubled to account for neglected contributions 
to the total gas mass (lower right panel). If a conversion fac- 
tor of Q = 4.6 (K km s~i pc'^)~i is used, both models can 
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Figure 3. N represents the number of halos one expects to find with mass ^ ^b,min/^b at a redshift z ^ 5 within a soUd angle 
of ^HDF = 4.4 arcmin^; for Af < 1, 1-N is thus the confidence with which an LVG model can be ruled out as a function of the 
minimum baryonic mass derived from the model (solid curve). Upper left panel: The vertical lines represent the mass values obtained 
for each of the four models presented in §3: (a) separate and unvirialized (dashed), (b) separate and virialized (dash-dot), (c) mixed and 
unvirialized (solid), and (d) mixed and virialized regions (dotted). In models (a) and (b), Mi,^j„ = Mj:/^ while in models (c) and (d), 
min = M/fj -|-M^ . Upper right panel: The minimum baryonic masses obtained for each model were doubled to account for neglected 
contributions to the total gas mass; models (b) and (d) can now be ruled out at the 2.5(T and 2.8(t levels respectively. Lower left panel: 
The vertical lines represent the mass values implied by the predicted CO(l-O) line luminosities from models (a) and (b) with a CO-to-H2 
conversion factor of o = 0.8 and 4.6 Mq (K km s"'^ pc^)"'^. If these minimum baryonic mass values are then doubled {lower right panel), 
both models can be ruled out at the > 2a level in the case where a = 4.6 Mq (K km s~^ pc-^)~^ is adopted as the conversion factor. 



be ruled out at the ~ la level and increasing these obtained 
masses by a factor of two drives up the sigma levels to > 2a 
for both models. 



5 SUMMARY 

In this paper, we employed the LVG method to explore alter- 
nate model configurations for the CO and C"*" emission lines 
regions in the high-redshift source HDF 850.1. In particular, 
we considered emissions originating from (i) separate unviri- 



alized regions, ('iij separate virialized regions, ('m^ uniformly 
mixed unvirialized regions, and (iv) uniformly mixed virial- 
ized regions. The set of LVG parameters, {Tkin, nH2, dv/dr, 
Xco' ^^2}, were fit to reproduce the observed line ratios 
and the resulting column density was then used to calculate 
the molecular gas mass, , in each case. In the models 
where the CO molecules and C"'' ions were assumed to be 
uniformly mixed, two additional parameters, nn and Nh, 
were fit as well, resulting in both a molecular and an atomic 
hydrogen gas mass. The gas masses derived by employing 
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the LVG technique thus represent conservative estimates of 
the minimum baryonic mass associated with HDF 850.1. 

These estimates were then used, together with the 
Sheth-Tormen mass function for dark matter halos to cal- 
culate the number of halos with mass Mh ^ Mb,,nin that 
each model predicts to find within the HDF survey volume. 
Given that at least one such source has been detected, we 
found that models (ii), (in), and (iv) can be ruled out at 
the 1.6a, 0.5a, and 2. Oct levels respectively. The confidence 
with which these models are ruled out is significantly in- 
creased if a less conservative estimate of the baryonic mass 
is taken; increasing the LVG-modeled gas masses by a fac- 
tor of two to account for neglected contributions to the to- 
tal baryonic mass, drives up these sigma levels to 2.4ct, 1.6ct, 
and 2.8ct respectively. We are therefore led to the conclusion 
that HDF 850.1 is modeled best by a collection of unvirial- 
ized molecular clouds with distinct CO and C"*" layers, as in 
PDR models. The LVG calculations for this model yield a 
kinetic temperature of 85 K, a velocity gradient of 5.0 km 
s~^ pc~^, a molecular hydrogen density of lO^ '^ cm~^, a CO 
to H2 abundance ratio of 10~^'*, and a column density of 
5.1x10^* cm~^. The corresponding molecular gas mass ob- 
tained using this LVG approach is Mj^^ « 2.54x10^^° Mq. 
For this preferred model we find that the CO-to-H2 lumi- 
nosity to mass ratio is a = 0.6 (K km pc'^)"^, close to 
the value found for ULIRGs in the local universe. 
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